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We consider low rank density operators q supported on a 
M x N Hilbert space for arbitrary M and N (M < N) and 
with a positive partial transpose (PPT) q Ta > 0. For rank 
r (o) < N we prove that having a PPT is necessary and suf- 
ficient for q to be separable; in this case we also provide its 
minimal decomposition in terms of pure product states. It 
follows from this result that there is no rank 3 bound entan- 
gled states having a PPT. We also present a necessary and 
sufficient condition for the separability of generic density ma- 
trices for which the sum of the ranks of g and q Ta satisfies 
r(g) + r(g TA ) < 2MN - M - N + 2. This separability con- 
dition has the form of a constructive check, providing thus 
also a pure product state decomposition for separable states, 
and it works in those cases where a system of couple polyno- 
mial equations has a finite number of solutions, as expected 
in most cases. 
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I. INTRODUCTION 

Entanglement is one of the quantum properties with 
no classical counterpart. It is closely connected to fun- 
damental questions of quantum mechanics an< ^ to 
physical phenomena which are important for quantum 
information processing [gj. The relevance of entangle- 
ment effects was first demonstrated for pure states. How- 
ever, in realistic physical situations one deals usually with 
mixed states, in which pure state entanglement has been 
significantly weakened by noise. In order to overcome 
the problems caused by noise (i.e. in order to reduce it) 
the idea of distillation of entanglement in spatially sepa- 
rated laboratories was introduced Q . It was proved then 
H that for bipartite systems of low dimensional Hilbert 
space C AI x C N or simply M x N (namely systems with 
M = 2 and N = 2 or 3) mixed state entanglement can 
always be distilled into its pure form. However, it turned 
out that in higher dimensional systems (MN > 6) bound 
entanglement |g] - which cannot be distilled, as opposed 
to free entanglement - exists. 

Unlike in the case of pure states, it is in general very 
difficult to know whether a given mixed state is entangled 
(inseparable) or non-entangled (separable). According 
to the definition a state supported on a Hilbert space 
T~Cab = Ti-A <8 Hb is separable if and only if it can be 
written in (or approximated by) the form Q 



Q = ^2pi\e l ,f l )(e l ,f l \, ^2pi = l, (1) 

i— 1 i 

where \&i,fi) stands here for a normalized vector je^) ® 
\fi) G Ha 8> Hb- In finite dimensional cases, the ones 
we will be concerned here, the approximation part is not 
necessary, as for any separable state one can always find a 
set {\ei, fi)} of product vectors for which k < dim(7i J 4s) 2 
in the above formula 

Several necessary conditions for separability are 
known: Werner derived a condition based on the analysis 
of local hidden variables models and mean value of the, 
so called, flipping operator 0; the Horodeckis proposed 
a necessary criterion based on a-entropy inequalities |j) . 
Peres demonstrated that the partial transpose q Ta of the 
matrix g, defined as (to, h\q Ta \n, v) = (n, ^|p|to, v) for 
any fixed product basis | n,ij) — 16^)^4 ® \&v) b > must be 
still a legitimate density matrix if g is separable Jl0| . 
This operationally friendly, necessary condition, called 
positive partial transpose (PPT) condition, turned out 
to be very strong. 

Soon after Peres result, a general connection between 
positive map theory and separability was established in 
pd| , where necessary and sufficient separability condi- 
tions were derived in terms of positive maps. In par- 
ticular it implied that for systems of low dimensions 
(MN < 6) the PPT condition is also sufficient for sep- 
arability. It was also shown that this is not the case 
for systems of higher dimensions (MN > 6). Later on 
explicit counterexamples of entangled states with PPT 
property were provided by means of another separability 
criterion, based on the analysis of the range of the density 
matrix || (cf. Jl^]). It was then shown that they repre- 
sent bound entanglement Q . Let us note that on math- 
ematical grounds there were examples, provided earlier 
|]l3| , of elements of positive matrices cones which can be 
treated as prototypes of PPT entangled states. 

Sufficient conditions for separability are also known. 
We remark that the results of |l4| readily imply that 
any state close enough to the completely random state 
7r is separable. Thus, as quantified in jj^ j, any mixture 
g = (1 —p)g + pit in a M X N system is separable if 
p > (1 + 2/MN)^ 1 or, in other words, as we wish to 
make explicit here, a full rank mixed state is separable 
provided its smallest eigenvalue be greater than or equal 
to (2 + MiV)- 1 . 

On the other hand the analysis of the range of the den- 
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sity matrices, first applied in the separability criterion 
f|, led to an algorithm for the optimal decomposition 
of mixed states into a separable and an inseparable part 
p6| , and to a systematic method of constructing exam- 
ples of PPT entangled states and peculiar positive maps 
0,^8|. Also, the technique of diminishing the rank of 
a PPT density matrix by subtraction of selected prod- 
uct vectors, which was worked out in [jl9), turned out to 
be very useful. This and other techniques have allowed 
quite recently to study operational necessary and suffi- 
cient separability conditions for states of a 2 x N system 
. In particular it has been shown that: 
(i) all PPT states of rank smaller than N are separable; 
(it) the separability of generic states such that r(g) + 
f(0 TA ) < 3./V reduces to analyzing the roots of some 
complex polynomials (a constructive separability crite- 
rion was derived, thus providing also the decomposition 
of such separable states into pure product states); 

(Hi) states invariant under partial transpose, and those 
that are not "very different" from their partial transpose 
are necessarily separable. 

This paper can be considered an extension and gener- 
alization of Ref. pp[ . The results (i) and (ii) obtained 
there for 2 x N systems are here generalized non-trivially 
to the case of M x N systems (M < N). We show, 
namely, that 

• any state g supported on MxN (M < N) and with 
rank r(g) < N is separable iff its partial transpose 
is positive; 

• separability of generic PPT density matrices with 
r(g) + r(g TA ) < 2MN - M - N + 2 reduces to 
solving a system of coupled polynomial equations. 

In both cases a pure product state decomposition for 
separable states is obtained. 

Throughout this paper we make use of the following 
definition: we say that a state p acting on M x N is sup- 
ported on M x N if this is the smallest product Hilbcrt 
space on which p can act. Let us introduce the local ranks 
t(qa) and r(gs), where gA,B = T^a,bQ are the reduced 
density operators. It immediately follows from the first 
of the above results that there is no PPT bound entan- 
glement of rank 3. Indeed, a rank 3 state p either has 
at least one of the local ranks t(pa) and r(ps) greater 
than 3, and in this case is distillable |2q] (i.e., p TA is not 
positive) , or else can be supported on a MN < 6 or on a 
3x3 system, and thus is separable. This implies in par- 
ticular that the bound entangled states constructed with 
the UPB method and those based on the chess-board 
structure of eigenvectors pl[ | are optimal with respect to 
their ranks. 

For our second main result, concerned with those PPT 
density matrices for which the sum of ranks satisfies 
r(g) + r(g TA ) < 2MN — M — N + 2, we identify the 
eligible product vectors (that is, those that can appear 



in decomposition ([j]) if p is separable) with the solutions 
of a system of coupled polynomial equations. We an- 
alyze these equations, which are arguably expected to 
have only a finite number of solutions. For this case we 
present a constructive (i.e. leading to a product state de- 
composition) method to check separability. Also for the 
same case we discuss an alternative, constructive method 
to check separability numerically. These checks represent 
a necessary and sufficient condition for separability. 

We wish to remark the importance of having separa- 
bility conditions for low rank density matrices, especially 
in relation to unsolved problems concerning the nature 
of bound entanglement (BE). Notice that such condi- 
tions are of great value when trying to construct states 
with BE. Among the open questions we encounter the 
existence of BE having a non-positive partial transpose 
NPT (sec 0]). Also, whether a finite or a vanishing 
amount of free entanglement is required to asymptoti- 
cally create bound entangled states. There are in addi- 
tion several conjectures concerning bound entanglement 
(see [p| Jl^j2^ -p5||) among them the ones connected to ca- 
pacities of quantum channels and bound entanglement 
assisted distillation. Finally, we have been recently able 
to establish a general connection between low rank bound 
entangled states and positive maps. This connection al- 
lows for a systematic construction of independent linear 
maps in arbitrary dimensions, including 2 x N, where 
the procedures based on unextendible product bases do 
not work [H . The discussion of this connection will be 
presented elsewhere. 

This paper is organized as follows: we start by gener- 
alizing some needed results of |2(| related to diminishing 
the rank of p by subtracting projectors on product vec- 
tors; in Section III we present our theorem about the sep- 
arability of states with rank < N; in Section IV the nec- 
essary and sufficient separability conditions for generic 
matrices with r(g) + r(g TA ) < 2MN - M - N + 2 are 
formulated, and discussed in the context of 3 x 3 systems; 
finally, Section V contains our conclusions and acknowl- 
edgments. 

II. DIMINISHING OF THE RANK - 
GENERALIZATIONS 

Before we turn to the main results of this paper we 
need to generalize some of those presented in [^0| . 

Consider a state g of a M x N system satisfying 
q t a > o. Throughout this paper K(X), R(X),k(X), and 
r(X) denote the kernel, the range, the dimension of the 
kernel, and the rank of the operator X, respectively. By 
{\ a i))iLi and {\bi)}iLi we will denote orthonormal basis 
in Ha and TLb, and by |e*} we will denote the com- 
plex conjugated vector of |e) in the orthonormal basis 
1 1} a j • ' ' j \M) a in which we perform the partial transpo- 
sition; that is, if |e) = J2i=i then |e*) = J2i=i a t N)- 
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For the time being we do not require M < N. The 
following Lemma is a generalization of Lemma 6 of Ref. 
{H proved there for M = 2: 

Lemma 1.- If 3 \ f) £ C w such that |a i; /) € for 
i = 1, M — 1, then either (i) |om, /} G -K"(f?) or (ii) 



Q\dM,f) = \a.M,g) 
TA \aM*,f) = \a M *,g) 



(2) 



for some \g) £ C N . 

Proof.- From the assumptions we have immediately 
Q TA \ai*,f) = (i = l,...,M - 1). In particular 
V|/i) £ C N we have (au*, h\g TA \a,i* , f) = or, equiv- 
alently, (a,-, h\g\a,M, f) — 0. Since \h) is arbitrary we 
have either statement (i) or g\a,M,f) = \o-M,g) for some 
\g) ^ 0. The second case needs further analysis. In 
a similar way we can prove that either q Ta \clm* , f) = 
(which is still equivalent to the statement (i)) or 
q Ta \o,m* , f) = \a>M*,g') for some \g') ^ 0. It remains to 
prove that \g') = \g). Indeed \g') = (a M *\Q TA \a M * , f) = 
(a M \g\a M i f) = Iff)- 

The second lemma below is also a generalization of the 
results from Ref. @: 

Lemma 2.- If g satisfies the assumptions of Lemma 
1, and the possibility (ii) of Lemma 1 holds, then 



Q\ = q- M a M,g)(aM,g\, 

where A = (a M , glg^ 1 \a>M, g) and 



(3) 



(i) g\ is a PPT state with r(gi) = r(g) — 1 and r(g^ A ) — 
r(g TA ) - 1. 

(ii) Qi is supported either on a (M— 1) x (N — 1) or on a 
M x (N-l). 

(iii) g\ is separable iff g is separable. 

Proof .- Following Corollary 1 and Lemma 2 from Ref. 
|E0| , we observe that 



qi = g - 



\aM,g)(aM,g\ 
(0M,5|e -1 |aM,fl') 



(4) 



is positive, and that r(gi) — r{g) — 1. Then (i) follows 
from taking into account that since 



A" 4 = {a M ,g\g 1 \a M ,g) = 
(g\f) = {a>M*,g\Q TA \a-M*,g}, 

we have that 



Qi —9 



\aM*,g)(aM*,g\ 
(aM*,g\e TA ~ 1 \a M *,g) 



(5) 



(6) 



so that also gJ A is positive andr(p^ A ) = r{g TA ) — 1. From 
assumptions on the kernel of g it follows that all vectors 



{|a i ,/)}f£ 1 belong to the kernel of qi, hence gi can be 
embedded into a M x (iV — 1) space; on the other hand 
since gi^A = QA — Mg\g)\ a M)(dM\, t {qi,a) must be either 
M = t{qa) or M — 1, which finishes the proof of (ii). 
In order to prove (iii) let us assume that g is separable 
and let us show that also Q\ is so (if g\ is separable then 
obviously g is also separable). Since g\a i} f) = (i = 
1, M — 1), we can always write 



(7) 



where (f\fi) = and 77 is a positive operator acting on 
C . If we impose \a,M,g) = Q\ a M,f) we obtain \g) = 
i]\f), and therefore |g) G R(rj)- We can write 

Pi = \ e i>fi)( e i>fi\ + l a A/)( a M| ® (l) - A|.g)(c/|), (8) 

so that if we show that the operator (rj — X\g}(g\) > 
then we have that p\ is separable. Using (g) we have 
that such an operator is 

v -m m=v -jgW^ l9){91 (9) 

and that therefore it is positive (cf. Lemma 1 in p0|). □ 



III. ALL RANK N PPT STATES SUPPORTED ON 
A M x N SYSTEM (M < N) ARE SEPARABLE 

In this section we generalize the following theorem 
proved in Ref. (§0|: 

Theorem.- [Theorem 1 of Ref. @] 

Let g be a PPT state of rank N supported on a 2 x N 
space. Then g is separable and can be written as 



N 



= \ e i>fi)( e i'fi\ 



(10) 



1=1 



with all {\fi)} linearly independent. 

We will express a density matrix in terms of its reduced 
operators (m|£?|j'a) acting on Hb- For instance, we will 
write (|l(]) as 



A 2?t 
B C 



(11) 



where A = (l\p\l) > 0, C = (2\p\2) > and B = (2\p\l). 
More generally, a M x TV density matrix will be expressed 
as 



En E12 

E\2 ^22 



E 



E 



IM 



where now E», 



(uleli/ 

ous Theorem to prove: 



ltd 

Z Z ' (12) 

... Emm _ 

We start by using the previ- 
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Lemma 3.- Let g be a rank TV PPT state supported 
on a 2 x TV space. Then after a reversible local filtering 
operation [] the state is proportional to the matrix (called 
hereafter 2 x TV canonical form) : 



B^B St 
B I 



[B I?[B I], 



(13) 



with B normal, i.e. [B, B'] = 0. 

Proof .- We write the density matrix (jl^) in the form 
(pi]). Because there is only a finite number of |e») in jlC| ) 
we can always find a vector \a) such that (a\ei) ^ for all 
i. Let this \a) be the second element of the orthonormal 
basis in Alice's space, i.e. 2} = \a). The matrix 

N 

C = (a\g\a) =Y J M^)\ 2 \U)(U\ 



has then maximal rank TV, since the \fi) are linearly in- 
dependent. Taking the local filter V — (VC*)" 1 on Bob's 
side (this corresponds to sandwitching the state between 
I (8> V and / ® V') we obtain 



Q 



A Bt 
B I 



(14) 



which is still positive and PPT (because any local oper- 
ation preserves the PPT property M). We can write 



e = E + diag[A,0], 



(15) 



where the positive matrix £ from the expression @) has 
rank r(£) = TV as its kernel K(£) has at least dimension 
TV containing all vectors of the type 



1^ = 11)1/) + 12)1-5/), 



(16) 



while its range has at least dimension TV due to the 
identity entry on the diagonal. Notice that diag[A, 0] 
is also positive, because positivity of g implies that A = 
A - B^B > (§7). Now, since in addition r(p) = r(£), 
we also have that R{p) = R(E) 2 i?(diag[A, 0]), that 
is if(diag[A,0]) 2 K(E). But K(E) is spanned by the 
states (pif), for which then (</>/|diag[A,0]|</>/) = 0, which 
finally implies A|/) = V|/). This ends the proof of 
the fact that A = 0, or in another words that A = B^B. 
This proves therefore the canonical form (|l^), but not 
yet the normality of B. The latter property can be sim- 
ply proven from the positivity of q Ta , which implies that 
RBt - B^B > p7|. The latter (positive) operator has 



1 A reversible transformation is a transformation that can be 
reversed with nonzero probability. A local filtering in Bob's 
side I®V is then reversible iff the operator V can be inverted, 
i.e. iff V' 1 exists. 



at the same time null trace and therefore it must vanish. 
Thus B is normal as stated. □ 

Let us prove now the generalization of Lemma 3 to the 
case of 3 x N systems (N > 3), and then to the M x N 
case, where M < N from now on. 

Lemma 4.- Let g be a PPT state of rank N in a 
3xJV space. Let the reduced state gs and the entry £33 
in some local basis have also the same rank TV. Then g 
can be transformed using some reversible local transfor- 
mation to the canonical form: 



q~\C,B,1$\C,B,I\ 



(17) 



where C, B are normal and [B, C'] = [B,C] = 0. 

Note that in Lemma 4, in contrast to Lemma 3, we 
assume that in some basis r{E^) = N. Later on, in 
Theorem 1, we will prove that this assumption is always 
satisfied. 

Proof . - In order to obtain the identity matrix / at the 
diagonal we use an analogous reversible local filter to the 
one used in the proof of Lemma 3. After that we readily 
obtain the form 



D B^B fit 
C B I 



(18) 



with both B and C normal and some unknown D. In- 
deed, expression ( |l8|) as well as the normality of B and 
C follow from the fact that after a local projection by 
projectors P k ®I= {\k)(k\ + |3)(3|) ® J, k = 1,2 we get 
a 2 x TV state satisfying the assumptions of Lemma 3. 

Now, notice that (H> f \g\y f ) = for |*/) = |2)|/) - 
\3)\Bf). Since g > we have that 



\l)\Dtf-&Bf), 



(19) 



which, as / is arbitrary, leads to D< — C^B. Thus for- 
mula (0) holds. Finally e shall use the latter as well as 
normality of B and C to prove that [B, C f ] = [B, C] = 0. 
We have 



0C B^C C 
C'^B B^B B 



(20) 



and we can check that for any |/) £ C N and for |<E>/) = 
|2}|/) - |3)|Bt/), {<S> f \Q T *\<f> f ) =0. As g is PPT this 
implies that 



g T ^ f ) = \l)\[B\C]f) 



(21) 



must vanish. Since the above equation holds for arbitrary 
I/} we have immediately that [B,&] = [C,B^ = 
Normality of B and of implies that these operators 
can be decomposed as a complex linear combination of 
projectors into eigenvectors. That they commute means 
that they actually have the same eigenvectors, and thus 
so do B and C, i.e. [B, C] = 0.D 
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Lemma 5.- Any PPT state supported on a M x N 
space (M < N) satisfying that (i) r(g) — N, (ii) in some 
product basis r{En) = N for some i, can be transformed 
after a reversible local transformation to the canonical 
form: 



Q 1 



Z^Z =[C 1 ,...,C M -iJ} i [C 1 ,...,CM-i,I} (22) 



with [C^C]] = [d, Cj] = 0, i,j = 1, M - 1. 

Proof .- It follows easily from the application of Lem- 
mas 3 and 4. In particular one has to use the local pro- 
jections P ® I = {\k)(k\ + \M)(M\) <g) I, 1 < k < M, 
P> (g> I = (\m){m\ + \m')(m'\ + \M){M\) ® I, 1 < m < 
m! < M.a 

As an immediate consequence we have 

Lemma 6 .- Any PPT state supported on a M X N 
space (M < N) satisfying that (i) r(g) = N, (ii) in some 
product basis r(Eu) = N for some i, is separable and 
can be expressed as 



JY 



^ fi){ e iifi\i 



(23) 



where the {|ei}} are possibly unnormalized and the 
are linearly independent. 

Proof .- We make use of Lemma 5. It is easy to see 
that the matrix Z^Z has nonzero eigenvectors of the form 
|oj,6j). Here \bi) is the i— th common eigenvector of all 
operators Cj, Cj, while \ei)^ — [cj , c^ -1 ^, 1] is a 
row of all i— th eigenvalues of matrices C\, Cm-i, I- 
Thus, after some reversible local transformation the state 
p becomes: 



N 



Q - 



^2 \ e i' b i)( e i> b i 



(24) 



where the are orthonormal. Reversing the previous 

local filtering we obtain (|23} ) □. 

Remark .- The above procedure gives a constructive 
algorithm to decompose any state which satisfies the as- 
sumptions of the Lemma. 

The main disadvantage of the above results is that all 
of them contain assumptions about r(E^) = N for some 
i and for some product basis, which as we have men- 
tioned, are not necessary. Our main theorem is free of 
that assumption (i.e., it shows that such En always ex- 
ists). To prove it we have to use induction with respect 
to M + N = K and use the previous Lemmas. We con- 
sider only r(g) = N, as a PPT state supported on M x N 
cannot have smaller rank. Indeed, since r^gs) = N, if 
r(g) < N then g is distillable, which implies that g TA is 
not positive pq ]. 

Theorem 1.- All rank N PPT states g supported on 
M x N are separable. 

Proof - We will prove that in some product basis we 
have r(En) = N for some i. Separability of g will fol- 
low from the previous Lemmas. Let us observe that 



the Theorem and the latter fact are true for M = 2 
and arbitrary N > 2. In particular they are true for 
M + N = K = 4 and 5. Let us assume that they hold for 
M + N < K, We shall now demonstrate that they also 
hold for M + N = K + 1. 

To this aim let us consider the case of g supported on 
a (M + 1) x N space, with M + 1 < N, r(g A ) = M + 1, 
t{qb) — N. and M + N = K. In an orthonormal, 
product basis representation the state g has the form of 
a (M + 1) x (M + 1) matrix with N x N entries 



Q = 



En 

eU 



El2 

E22 



Ei, 



M+l 



E 



l.M+l 



Em+\,m+\ 



(25) 



Let us consider the following M x M submatrix of g 



W(g) = W 



E 2 
E 



■2:; 



E23 
E22 



E- 



2.M+1 



E. 



t 

2, M+l 



E 



M+l, M+l 



(26) 



resulting after removing the first row and the first column 
from the representation ( p5[ ) . As the latter action can be 
achieved by a local projection on Alice's side, W is an 
unnormalized PPT state acting in M x N . For Bob's 
reduced matrix Wb = TrB(W), we shall consider two 
alternative possibilities: (i) t(Wb) = N, (ii) r{Ws) < N. 

In case (i) we must have r(W) = N as otherwise we 
would have that the global rank is less than one of the 
local ranks, resulting in distillability of g, ergo in viola- 
tion of the PPT condition |],^(|. But it means that W is 
a PPT state supported onMxf with global and local 
rank equal to N. According to the induction assumption, 
it is thus separable and for some product basis has an en- 
try En for some i = 2, . . . , M + 1 with the rank N. But 
then g has an entry En with rank N in the same product 
basis, and from Lemmas 5 and 6 it follows immediately 
that it is separable. 

Consider now case (ii). If W has t{Wb) < N, then ob- 
viously there exist a sequence of product vectors \at, f) £ 
K(W), i = 2, ...,M + 1. It is immediate to check that 
they must belong to kernel of g. That means that the 
assumptions of the Lemma 1 are fulfilled. The possibil- 
ity (i) of this Lemma cannot hold because otherwise one 
could embed g into [M + 1) x (N — 1) space, and r{gs) 
would be N — 1 instead of N. The possibility (ii) of 
Lemma 1 means that g can be written in the form (cf. 
Lemma 2) 



Q = 3 



A|l,<?)(l,5l 



(27) 



where g 1 is a rank N — 1 PPT state supported either 
on a (M + 1) x (N — 1) subspace or on a M x (TV — 
1) subspace, A -1 = {l,g\p~ 1 \l,g) and \l,g){l,g\ is an 
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unnormalized projector onto a product state such that 
1 1 1, g) is orthogonal to R(g'). 

At the same time it must hold that r{g' B ) = N — 1, 
since i) Bob's space has now only N — 1 dimensions, ii) 
r(g' B ) cannot be smaller than N — 1, since N = r{gs = 
Q'b + IfXfl) an d IffXffl can increase at most in one unit 
the rank of g' B . All that means that the matrix g' fulfills 
the induction assumption as (M + 1) + (N — 1) = K 
(or M + (N - 1) = K - 1 ) and r(g') = r{g' B ), ergo it is 
separable and has in some product basis a^, bj) the entry 
E' u = (a,i\g'\ai) with rank N — 1. Lemma 6 implies then 
that g (= g' + \\l,g)(l,g\) can be decomposed into 



E 

2=1 



e 4 , fi)(ei, fi\ + X\l,g)(l,g\, 



(28) 



where \g) is linearly independent from the set of (also 
linearly independent) vectors Since there is only a 
finite number of projectors in the decomposition above, 
we can always find a vector |o) in Alice's space such that 
{ei\a) ^ 7^ (l|e). Including such a vector in a product 
basis to express g we will obtain the wished rank N ele- 
ment (a\g\a). This completes the proof of the induction 
step, and by induction completes thus the proof of the 
theorem. □ 



IV. SEPARABILITY CRITERIA FOR 

RANK(5)+RANK(/") < 2MN -M-N + 2 

In this section we generalize the results obtained for 
2 x N systems in Ref. The idea is that a PPT den- 
sity operator g with r(g) + r(g TA ) < 2MN - M - N + 2 
may have a finite number of product vectors e^, fi) in 
its range, such that |e|,/j) G R{g TA ). These product 
vectors are the only possible candidates to appear in de- 
composition (|l|) Q . Finding them requires solving a sys- 
tem of polynomial equations. First we show how to solve 
these equations in a generic case (namely when the coef- 
ficients of such equations do not happen to satisfy a large 
scries of conditions, which amounts to having only a fi- 
nite number of solutions), and once all the product states 
{l e ii fi)}f=o° have been obtained, we present an algorith- 
mic method to check whether p is separable. This is done 
in a finite number of computational steps, and thus solves 
operationally the problem of separability for states with 
r(g) + r(g TA ) < 2MN - M - N + 2 and finite L. 



A. Eligible product vectors. 

Let the linearly independent vectors \Ki), \Ki) form a 
basis in the kernel of g and in the kernel of g TA , respec- 
tively: 



K(g) =span{|Jf i ),t= l,...,k(g)} (29) 
K{g TA ) = span{|X i ),i= 1, k(g TA )} (30) 

We consider here the case when k(g) + k(g TA ) > M + 
N — 2. We can always expand \Ki) and \Ki) in an or- 
thonormal basis in Alice's space 



\Ki) 



M 

E 

m=0 



m,k™), 



M 

\Ki) = J2 K^™}- 

rn— 



(31) 
(32) 



A product vector |e, /) belonging to the range R(g) must 
be orthogonal to all \Ki); simultaneously, if its partial 
complex conjugation belongs to R(g TA ), \e* , /} must be 
orthogonal to all \Ki). Thus the eligible product vectors 
are the solutions of k{g) + k(g TA ) equations, namely 



(Ki\e,f) = 0,i = l,...,k(g), 
(Ki\e*,f)=0, i=l,...,k(g TA ). 

Let us now expand e) in the above formula as: 



|e) 



Oil 



Um 



(33) 



(34) 



We restrict ourselves to a\ = 1. The reason is that 
we expect to find only a finite number L of inequivalent 
vectors \ei, fi) that fulfill the requirements. A generic 
choice of an orthonormal basis { |a^) } in Alice's space will 
imply that (l\ai) ^ for all i = 1, . . . , L. In this basis 
«i can be set equal to 1. 

Equations ([33]) can be rewritten as follows: 



A(ai, ...,a M ;a*i, -,a* M )\f) = 0, 



(35) 



where the [k(g) + k(g TB )] x N matrix A is defined as 
follows: 



ai, . 



Ef=i"™( fc fc 



(e) 



Ei=l a m(k^i T A 



k{ S TA V 



D k ( e )y N {a>) 



D 



k( e T B)xN 



(«•) 



(36) 



If ( p5| ) holds for some \f) ^ and |e) ^0, this means 
that for the corresponding set of a's the rank of A is 
smaller than TV. Therefore, in order to identify eligi- 
ble product vectors we have to require that at most 
N — 1 rows of A be linearly independent vectors. In 
what follows we restrict ourselves to the limiting case 



G 



k(g) + k(g TA ) — M + N — 2, the others containing more 
restrictions and consequently less solutions than this. 

Let us then take N — l rows of A, say the first ones, and 
let us require that each of the remaining M — 1 rows be 
linearly dependent of them. Recall that we can use the 
M— 1 variables aa, • • • , olm in order to achieve this. Then, 
parameter counting strongly suggests that we need to fix 
all the M — 1 a's in order to make A have rank smaller 
than N, this corresponding to a zero measure set of points 
in the a-space [ai = l, ■ ■ ■ , ckm]- We will in addition 
relate the number of solutions to the roots of complex 
polynomials, which under generic conditions have only 
a finite number of roots. Numerical experience acquired 
for the 2 x N case further supports the expectation that 
the number of solutions be typically finite. 



B. Generic Polynomials. 

Let us discuss a bit further sufficient conditions for the 
existence of a finite set of solutions, while presenting a 
systematic method to find them once the conditions are 
fulfilled. This method also works for k(g) + k(g TA ) < 
M + N — 2 by just adding more equations. 

Matrix A will have at most rank N — l after requiring 
that all its rank N minors vanish. At risk of finally find- 
ing more solutions than just those of equations (35), we 
can impose that only M — 1 of these minors vanish. The 
reason for doing so is that this will already allow us to 
prove that only a finite number of product vectors fulfill 
( |33| ) under some generic circumstances. Thus we con- 
sider the determinant of N x N submatrices of A formed 
by taking its first N — 1 rows and then also one of the 
M — 1 remaining ones. We shall denote these minors by 
Fi{z\, . . . , Z2m), i = 1, • • • ,M — 1, where Zj = aj and 
zj +M = a* (i — 1, . . . , M) will be taken as 2(M - 1) 
independent variables {z\ = zm+i = 1)- Again, this will 
only imply that when we now set 



FAz!, 



,Z 2 M 



0, 



(37) 



for i = 1, • • • , M— 1, some of the solutions we find do not 
correspond to product vectors, although all the \et, fi) 
we look for are among the solutions of (|3^) . 

We have 2{M — 1) variables Zi and the same number 
2(M — 1) of polynomial equations for them, M— 1 coming 
from the minors Fi(z\, Z2m) = and the remaining M — 1 
from its complex conjugation, which are inequivalent to 
the first ones as the variables are mapped according to 
Zi <-> z.i + M, i = 1) • • • > M, under complex conjugation. 

No theorem exists for complex polynomials P(a,a*) 
which allows us to know the number of roots they have. 
However, in a generic case, namely when P(a, a*) is not 
proportional to its complex conjugate, we can prove that 
only a finite number of solutions exist. In p(i| a method 
to find such roots was developed for polynomials depend- 
ing on one a and its complex conjugate. Accordingly, 



from P(a, a*) another polynomial Q(a) containing all 
the roots of P was obtained. Such a method admits a 
generalization to the present case, which we shall discuss 
later on by means of an example when analyzing states of 
a 3 x 3 system. As already mentioned, we were not able 
to determine when a density matrix g will lead to a set of 
non-generic polynomials. However, we expect this to be 
rarely the case. In what follows we will assume that the 
polynomials derived from g are generic, and that there- 
fore there is only a finite number of product vectors that 
can appear in (m). 



C. Separability criterion. 

When the number of solutions of equation ( ^3|) is finite, 
we can formulate a necessary and sufficient separability 
condition which follows from the following general theo- 
rem: 

Theorem 2 (see also ||J) .- A state g of rank r(g) 
is separable iff it can be written as a convex combina- 
tion of at most mm{r(g) 2 ,r(g TA ) 2 } linearly independent 
projectors |e», fi)(ei, fi\ onto product vectors. 

Proof.- The inverse implication is obvious. For the 
direct implication we will assume, without loss of gener- 
ality, that r(g) < r(g TA ). Caratheodoris' theorem tells 
us then that g can be expressed as a convex combination 
of r(g) 2 product projectors, 



r(ef 

Q = ^2 Pi\ei, fi){ei, fi\ 



(38) 



Suppose these projectors are not linearly independent. 
This means we can find J^i c i\ e ii fi){ e ii fi\ = with at 
least some non- vanishing Cj S 1Z. Set A = min{pi/ci}. 
Then the decomposition 



r(e) 2 



0= ^2 (Pi - *>Ci)\ei, fi)(ei, fi\, 



(39) 



corresponds also to a convex combination of the previous 
projectors | e^, /i)(ej, /,-|, but with at least one of the terms 
having vanishing weight. Now, if the remaining projec- 
tors do not form yet a linearly independent set, we can 
repeat the same procedure and get rid of another prod- 
uct projector. This can be iterated until expressing g 
as a convex combination of linearly independent product 
projectors. □ 

Consequently, once we obtain all product vectors 
\e h fi) € R(g) such that |e*, /,) e R{g TA ), i = 1, • • • , L < 
oo, we can find out whether g is separable by proceeding 
as follows: 

• We build all possible maximal subsets of linearly 
independent projectors \ei, fiYfii, fi\ (with at least 
Lq = max{r(p), r(g TA )} elements). Notice that 
there is only a finite number of subsets. 
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• For each of these subsets we express g as a linear 
combination of projectors in the subset. 

• If this is possible, then we have to see whether the 
coefficients of the linear combination are all posi- 
tive. 

We immediately have: 

Separability criterion.- g is separable iff all coef- 
ficients are non-negative in (at least) one of the linear 
combinations described above. 



D. Numerical methods. 

We notice that for a g with just a finite, but large 
number L of eligible product vectors it may be impracti- 
cal to construct all possible subsets of linearly indepen- 
dent product projectors, as described above. In this case 
the linear programming theory f2S[| has developed vari- 
ous methods to try to find out a solution to whether g 
can be expressed as a linear combination, with positive 
weights, of the over complete but finite set of projectors 
\ e ii fi)( e ii fi\- We propose, however, to use for this aim 
the best separable approximation (BSA) method, devel- 
oped by us in Ref. @. It has nice physical analogies 
also for non-separable states, providing the expansion 

g = g s + (1 - X)5g, 

where g s = J2 t AjP^ is a separable state, A = J2i * s 
maximal, and finally 5g is a state that does not have any 
product vector in its range. The paper p6| describes an 
efficient algorithm for finding such expansion, by opti- 
mizing each of the Aj individually, and each of the pairs 
Aj, Aj with respect to Aj + Aj. For the purpose of check- 
ing if a given matrix is separable, the BSA method of 
Ref. [[l6) is sufficient; in the context of the present paper 
it is interesting to introduce here a generalization of the 
results of @ to the PPT states ||: 

Lemma 7.- Let g be a PPT state. For a given set 
of P{ = fi)(ei, fi\, such that the product vectors 
\ei, fi) £ R(Sg), such that |e*, fi) e R(Sg TA ), there exists 
the best separable approximation of g, in the form 

Q = Qs + (1 - X)Sq, 

where g s = J2i AjPj is a separable state, A = Yli A» 
is maximal, and finally both Sg > 0, and Sg TA > 0. 
Moreover, there does not exist a product vector \e, /} 6 
R(5g), such that |e*,/} G R(Sg TA ). 

The proof of the above lemma is the same as the proof 
in Ref. |l6| . Similarly, one can find an efficient algorithm 
for finding the BSA, by requiring that: 

• All A, should be maximal, i.e. 



A, = min (ei,fi\(p -J2 A j P j) ^^iji) \ 



{e*,H(p TA - Y,kiPj A )- X \<Ji)- l \ ■ (40) 

All pairs of Aj,Aj should be maximized with re- 
spect to Aj + Aj. This requirement can also be 
expressed in an analytical form for A's, which will 
be presented elsewhere M. 



E. Example: 3x3 system. 

We end this section by describing with an example in 
a 3 x 3 system how to estimate the number L of eligible 
product vectors. This example illustrates how to gener- 
alize to several independent a's the method developed in 
& 

Suppose r(g) < 4 and r(g TA ) < 9. For r(g) = 4 and 
r(g TA ) = 9 (least favorable case) we have that the matrix 
A reads 



.4 



(fell 

(H\ 
(fell 



a 2 (k 2 1 \ 
a 2 {ki\ 
a 2 {kl\ 
a 2 {k\\ 
a 2 (kl\ 



a 3 (kf\ 
■a 3 {k%\ 

a 3 (kl\ 
■a 3 (kl\ 

a 3 (k!\ 



(41) 



so that after constructing the 3x3 submatrices ^1,2,3 by 
taking the first two rows of A and one of the remaining 
rows at a time, we obtain three 3-rd order equations for 
a\ and 



Fi = dct Mi = a 2 P 3 (a 3 ) = 0, 



k=0 
3 



F 2 = det M 2 = "2^3 (as) = 0, 



k=0 



F 3 = dct M 3 = 0, 



(42) 

(43) 
(44) 



where P s (x) denotes a s-th order polynomial in x. By 
only using equations ( ft^) and (^) we can obtain two 
quadratic equations in a 2 as follows: on the one hand 
we multiply @ by Q|(a 3 ), © by P 3 3 (a 3 ), and then 
subtract them, leading to 



k=0 



a k 2 R k 6 (a 3 ) = 0; 



(45) 



on the other hand we multiply @ by Q 3 (a 3 ), (H) by 
P§(a 3 ), again subtract them, and after dividing by a 2 
we obtain 
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2 

E 

k=0 



(46) 



Finally, applying the same trick but now to equations 
( [45| ) and (^6|), we obtain two linear equations for 0:2, 
from which a unique 18-th order equation for as fol- 
lows. Therefore there are at most 18 different values 
of 0:3 which in principle could lead to an eligible prod- 
uct vector. For each such values one should now still 
solve the 3 3rd order equations ([l2|-|44|) for 02, and see 
which solutions survive, if any p| . Finally, for those triads 
[1, a2, CX3] which indeed fulfill ( ff2|ff4| ) we can diagonalize 
A and find the corresponding Bob's local vector |/) in 
the kernel of A. We have obtained, thus, L < 18. 

Before going into the conclusions we shall discuss 
briefly the question of the relative size of r(g) and r(g TA ). 
It is natural to expect that this difference is not too big. 
However some naive intuitions must be abandoned (see 
pot). Here we shall make the simple observation: 

Observation .- Let g be a PPT state. If kernel of g 
contains the range of some PPT state a, then the kernel 
of g TA contains the range of <j Ta , so that r(g TA ) < MN — 
r(a TA ). 

The above observation about rank of g follows easily 
from the fact that Tr(AB) =Tt(A Ta B Ta ). Note that a 
can be a separable state. In particular, if the kernel of 
g contains any system of n orthogonal product vectors 
(in particular UPB set |Ti| ) then r(g TA ) can not exceed 
the value of MN — n. The same holds if a from our 
observation is PPT bound entangled state defined as a 
UPB complement |l7j. The rank of the latter does not 
change under partial transpose, so again r(g TA ) can not 
exceed the value of MN — r(<r). It can be also extended 
in other direction: taking a as a nontrivial PPT invari- 
ant state. Apart from all ct's being complements of real 
UPB's, there is an other nontrivial class (provided in |32|) 
of N x N states of that kind all having r(er) = +1. 
From the above discussion and the Theorem 1 we imme- 
diately know, for example, that for all the 3x3 PPT 
entangled states with the kernel containing UPB com- 
plement both ranks: r(g T2 ) and r(g T2 ) must amount to 
either 4 or 5 so they cannot differ much from each other. 



V. CONCLUSIONS 



rank. There are several problems, however, which remain 
open and are worth further studies: 

• In our analysis of the kernels of g and g TA we have 
essentially used only those of their properties that 
are consequences of the dimensionality. On the 
other hand, it is expectable that both kernels are 
structurally related through the partial transpose 
operation. It would be important to investigate 
such relations, since it would probably automati- 
cally put much more stringent restrictions on the 
existence of separable matrices of low rank. 

• All of the results of this paper can be generalized to 
the case of multipartite systems, and in particular 3 
partite systems. We have already obtained several 
results, but we leave a detailed and complete dis- 
cussion of this problem to a separate publication. 
Let us just mention here that according to our stud- 
ies: i) there are no rank N PPT entangled states for 
NxNxN systems; ii) In 2 x 2 x 2 spaces PPT states 
of rank 4 are separable with respect to 2 x 4 space 
of Alice and joint space of Bob and Charles, and 
posses a unique decomposition into a sum of 4 pro- 
jectors onto product vectors in 2 x 4 space; they are 
fully separable iff those 4 product vectors are at the 
same time product vectors in the sense of 2 x 2 x 2; 
iii) In 2 x 2 x 2 spaces generic PPT states with 
r(g)+r(g TA )+r(g TB )+r(g Tc ) < 4x8-2x2+1= 29 
have a finite number of product vectors in their 
range, such that the partial conjugates of those 
product vectors are in the corresponding ranges of 
partial transposes. 
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We have presented in this paper a relatively complete 
list of separability criteria for density matrices of low 



2 Notice that for a given 03 in principle we could find 0, 1, 2 
or 3 valid values of Q2. For simplicity we assume in the final 
estimation of the number L of eligible product vectors that 
to each solution Q3 there corresponds at most one valid 02. 
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